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r— 1 ^ Abstract. We argue, for a wide class of systems including graphene, that in the 

fH ■ low temperature, high density, large separation and strong screening limits the drag 

resistivity behaves as d^^, where d is the separation between the two layers. The results 

QJ ' are independent of the energy dispersion relation, the dependence on momentum of 

the transport time, and the electronic wave function structure. We discuss how a 
correct treatment of the electron-electron interactions in an inhomogeneous dielectric 
background changes the theoretical analysis of the experimental drag results of Kim 
et al (2011 Phys. Rev. B 161401). We find that a quantitative understanding of 
i-rt ' the available experimental data (Kim et al 2011 Phys. Rev. B 161401) for drag in 

^ ■ graphene is lacking. 
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1. Introduction 

Coulomb drag [H [2] occurs when driving a current through a metallic layer, referred 
to as the active layer (and which will be denoted by 2), induces a current in another 
metallic layer, separated by a distance d, referred to as the passive layer (and which 
will be denoted by 1). This phenomenon is caused by the transference of momentum 
between electrons in different layers due to the interlayer electron-electron interaction. 
In experimental situations no current is allowed to flow in the passive layer, so that an 
electrical field, Ei, builds up in that layer. In this situation, the drag resistivity, p^, is 
obtained from the ratio 

Pn = fi. (1) 

where J2 is the current driven through the active layer. 

Recently there has been great deal of interest on the phenomenon of Coulomb 
drag in graphene double layers. Although the number of experimental works on drag 
in graphene is still scarce |3l SI E], there are already plenty of theoretical works on 
the topic [El m El [9l [101 [HI Il2]- However, among the theoretical literature we find 
contradictory statements, particularly in the limit of low temperature, high density, 
large layer separation and strong screening. In this limit, it is a well established result 
that the drag resistivity between two 2DEG's (two dimensional electron gases), with 
parabolic dispersion relations and a constant transport time, depends on temperature 
as T^ and on the interlayer separation as d~^ [T3] (a result previously given in [H] 
without derivation). Graphene is distinct from the 2DEG in three ways: (i) the low 
energy dispersion relation is linear, instead of parabolic; (ii) the electron wave function 
has a spinorial structure, instead of being a scalar; (iii) the transport time for the 
dominant kind of impurities is proportional to the momentum, r(/c) oc k, [15]. There is 
agreement [3 El [9l [101 [HI US] that in the low temperature/high density limit, the T^ 
dependence should still hold in graphene. However, there are a few contradictory results 
on the dependence of the drag resistivity on the layer separation, and how considering a 
constant or a momentum dependent transport time changes the result. In this paper we 
will give special attention to the limit of low temperature, high density, large interlayer 
separation and strong screening. Therefore, we will make a brief overview of the available 
results in the literature in this limit. 

(i) Tse et al \l\ assumes a constant transport time, and a dependence of d~'^ is obtained. 

(ii) Peres et al [8] considers a momentum dependent transport time, T{k) oc fc, obtaining 
a d~^ dependence (as we explain later this result is due to an algebraic error at the 
end of the asymptotic calculation; correcting this gives a d~^ dependence). 

(iii) Katsnelson [9] considers a constant transport time and obtains a (i~^ dependence. 

(iv) Hwuang et al [10] considers both cases of a constant transport time and a transport 
time proportional to the momentum, T{k) oc k. For a constant transport time a 
d~^ dependence is obtained (in contradiction with the result from ^). For the 
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case of a momentum dependent transport time a d~'^ behaviour is obtained. The 
case of drag between two bilayer graphene layers is also studied. Using a constant 
transport time a d~'^ behaviour is obtained in this same limit. 

(v) Narozhny et al [H] considers both cases of a constant and a linearly momentum 
dependent transport time. For the case of a constant transport time it is obtained 
a c?""' dependence. It is argued that in the low temperature, high density limit, this 
result still holds, regardless of the momentum dependence of the transport time. 

(vi) Carrega et al [12], in a recent independent work, studies drag between massless 
Dirac electrons. It is proved that in the low temperature/high density limit, the 
dependence on momentum of the transport time is irrelevant and a d~'^ dependence 
is obtained for large interlayer separations. 

In this paper, we attempt to clarify this situation by presenting a clear proof that 
in the limit of low temperature, high density, large interlayer separation and strong 
screening the drag resistivity should always depend on temperature as T^ and on 
distance separation as d~'^. Our analysis is independent of the energy dispersion relation, 
electron wave function structure and dependence on momentum of the transport time. 
The structure of the paper is as follows: in section [21 we present the general theory 
of Coulomb drag. In section [3], we present a general argument proving that the drag 
resistivity in the limit of low temperature, high density, large separation between layers 
and strong screening should always depend on temperature as T^ and on the interlayer 
separation as d~'^. We also study the case of graphene for small interlayer separation. 
In section HI we specialize to graphene and re-analyze the asymptotic result derived in 
[8]. Finally, in section |5] we look into the experimental data of [3] using a more careful 
treatment of the bare Coulomb interactions. 

2. General formulation of Coulomb drag 

Theoretically, it is convenient to compute conductivities instead of resistivities. For 
isotropic systems the drag resistivity is related to the conductivities by 

PD = ^ , (2) 

0"llCr22 — Cri20"21 0"llO"22 

where an and (722 are the intralayer conductivities of the passive and active layers and 
a£) = cri2 = 0-21 is the drag conductivity; in (j2]) it was assumed that an, 0-22 ^ o"r>- 

Considering that tunnelling between layers does not occur and that the intralayer 
transport is dominated by impurity scattering, the drag conductivity can be computed 
in second order in the interlayer interaction using either Boltzmann's kinetic equation 
[TOl [H |13l [16], the memory-function formalism [T7] or Kubo's formula [TOl [HI [18], and 
is given by 

ij 6162 f d^q r+'^ du\Ui2{q,Uj)f ■ • 



On Coulomb drag in double layer systems 4 

where e^ is the charge of carriers in layer a, Ui2{q,u)) is the interlayer electron-electron 
interaction and r^(g', w) is the i-th component of non-linear susceptibility of layer a. In 
the weak impurity limit, the non-linear susceptibility reads [71 [TDj (the layer index is 
omitted for simplicity) 

72; 



-> V ■\ / d k „x y 



^F[%x)-^F[ej;^^y} 



k, X I k', X' 



is the 



X {^k,xTk,x - ^k+q,x'^k+,-,x') ^ [%x - ^k+<l,X' + ^ ) . (4) 

where g is the flavour degeneracy, A, A' are band indices, f-.'p 
electron wave function overlap factor (which encodes the structure of the wave function), 
riF^e) is the Fermi-Dirac distribution function, v^^^ is the particle's group velocity, tj:^ 
is the impurity transport time, and e^y^ is the energy dispersion. 

For the interlayer interaction one usually uses the RPA dynamically screened 
Coulomb interaction flGl. 



Ui2{q,uj) = -^— , (5) 

eRPA{q,uj) 

where Vab{q) is the bare Coulomb interaction between electrons in layer a and b, and 

^rpa{(1,'-^) is the RPA dielectric function for the double layer system, which is given by 

m 



eRPAiq,(^) = [1 - Vii{q)xiiq,uj)] [1 - V22iq)X2iq,uj)] 

-Vi2{q)V2i{q)Xi{q,(^)X2{q,Uj), (6) 

and Xa{q,(^) is the polarizability of layer a. The bare Coulomb interactions can in 



general be written as (see Appendix A ) 



Va,{q) = ^^e-'"^^^-^-), (7) 

eabiq) 2eoq 

where eo is the vacuum permittivity and eab{q) are effective dielectric functions. If 

the metallic layers are immersed in a homogeneous dielectric with constant e^ then 

eab{q) = er- 

3. Low temperature behaviour 

We now study the behaviour of the drag conductivity in the limit of low temperature, 
^Fi{2) 3> fc^T, where eFi(2} is the Fermi energy of layer 1(2). Unless specified otherwise, 
we will keep the energy dispersion relation, eg^, the transport time, Tj^ ^, and wave 
function overlap factors, /^.'^ _, general. We assume isotropy and that there is only 
one band at the Fermi level. Central to the analysis it the realization that the energy 
dispersion relation close to the Fermi energy is always linear in momentum, that is, the 
dispersion can be approximated by: 

ekc-ep^hvpik-kF) , (8) 
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where vp is the slope of the band at the Fermi energy, termed Fermi velocity, and the 
label c refers to the conduction band. For graphene ([8]) is exact. We also assume that 
the two metallic layers are placed in vacuum, such that 

(2 \ 2 
^— j 2 sinh (g(i) e"^"* 

-n — [Xi(9',w) + X2(g,w)]. (9) 

zeoq 

Due to the factor sinh~^ {(3hjj /2) in (|3]), the main contribution to the integral in 

u comes from tyjj < fc^T. Since ei7'i(2) ^ ksT, we can therefore expand the remaining 

integration kernel to lowest order in u and set T = 0. Therefore we replace the 

dynamically screened dielectric function, eijp^(g, w), by its static value, eijp^((f, 0), and 

expand the non-linear susceptibility of each layer, (|1]), to lowest order in uj. Using the 

energy conserving 5- function in (|1]), ^ \t^x^ ^k+qX' + ^j; we expand to lowest order 



m Lo: 



k,x + ^) 



^F [%X) - ^F [^k+<l,X') = ^^ {%X) - ^^ [^ 

^ -^ ^^ ^ ^^ [^F - 6g,,) . (10) 

Therefore T{q,u) has a linear contribution in uj. Since we want T{q,uj) to lowest order 
in UJ we can now set w = in 5 ( e^ ^ — eg_(_-.^, + ^ j , obtaining 

r(g, ^) = -91^ J d^kflij (ep - e^,)i 5 (e^ - 6,^^.^) 



X yn,Jk,c - n+q,Jk+q,c) ■ (n) 

Since we have isotropy we can write v^^t^^ = kg{k), where g{k) is a general function 
that satisfies kpg{kr ^ 



(5-functions set 



-F) 



vpTp, Tp being the transport time at the Fermi level. The 



k + q 



kp, and therefore we can take g{k) outside the integral, 
obtaining to lowest order in u 

Note that in this limit q is restricted to g < 2kp. To perform the integration in k, we 
choose, without loss of generality, q= (g, 0) and write 

- a ^'^ 
u = cost^ = — — , 

kq 

d^k=2kp r dk ^' '^'^ 



J-i a/1 -m^ 

6[ep-e,^^) = ^^6{k-kp), 

Siep-e^:^.) =-—6lu + 



^+'l'V hvpq V 2A; 
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Therefore the following result is obtained 



r(g,w) = 9 



Txhvp q 



rc,c 
•' k,k+q 



VT^ 



u 



(13) 



k=kp,u— 2_ 



2kp 

This result is central to this paper. It shows that in the limit of low temperature and 
high density, the non-linear susceptibility is independent of both the energy dispersion 
relation and the dependence of the transport time on momentum, depending only in 
the particular form of the overlap factor. Therefore, this result can be readily applied 
for the case of a 2DEG, graphene, bilayer graphene and other systems. Although it 
was already pointed out in [HI [T2] that the non-linear susceptibility is independent 
of the momentum dependence of the transport time, in those works this result was 
obtained for the particular case of massless Dirac electrons. Here, we show that this is 
a general result also independent of the energy dispersion relation. Note that this result 
contradicts [10] , where different results for the non-linear susceptibility are obtained for 
different transport times in the low temperature limit. 

Since in the low temperature limit we have r(g, u) oc cj, the integration in a; in ([3D 
reads 



/o sinh' (/3na;/2) V ^ / 6 ' 

which gives the T^ dependence of the drag conductivity and resistivity in the low 
temperature limit. The T^ behaviour is independent of the details of the energy 
dispersion relation, the transport time and the wave function overlap factors. Notice, 
however, that the T^ behaviour might be modified if one includes corrections to the 
drag conductivity due to higher order terms in the interlayer interaction [20] . 

3.1. General system at large interlayer distance and strong screening 

We now assume that the interlayer separation is large, kpd ^ 1. The interlayer Coulomb 
interaction decays exponentially with rf, thus the integration kernel of ([3]) is dominated 
by values of q such that q < d~^. Therefore the condition kpd 3> 1 allow us to expand 
the remaining integration kernel to lowest order in q. To lowest order, the overlap factor 
f^Z ^ is 1. Therefore at low temperature and for small q and cu, with u < vpq, we have 

rg,w =c/^ , (15) 

Trnvp q 

a universal result that is independent of all the details of the system. Note, that although 

it is clear that in this limit F should only depend on quantities defined at the Fermi 

level {kp, Tp), it is not obvious at first that changing the momentum dependence of 

e^ or r^ will not change the momentum dependence of F. For small q we approximate 

Xa{q,0) — —Pai^Fa), where Pa{^) is the density of states of layer a, and the RPA 

dielectric function ([9]) becomes 

/ n\ II lTFiqTF2^ ■ U / ^\ -qd , qTFl + qTF2 ,-.„. 

eijp^(g,0) = H 2smh(gd)e '^ H , (16) 
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with qxFa = Pa(eFa)e^/(2eo), the Thomas-Fermi screening momentum in 2D of layer a. 
If we assume that we have strong screening, qTFi{2)d ^ 1, we further approximate [13] 

enPAiq, 0) = 2^^^^ sinh (qd) e"'"^. (17) 



g2 

If the dispersion relation of layer a is given by a power law, e| = Cak^""-, then 

we have qrpa oc kp^ . Therefore, for a linear dispersion relation the condition 
qTFad >■ 1 is equivalent to kpad ^ 1; while for a parabolic dispersion relation qxpa is 
independent of kpa, and therefore qrpad ^ 1 becomes an extra assumption. Assuming 
QTFi{2)d 3> 1, and using (IT5I) and (TT7I) in Q we obtain the following expression for the 
drag conductivity: 

^ 6162 Ci^)9i92 ej el tfiTf2 {kpTf 

^ h 2'^ AiTeoVFihAireoVF2hn?{qTFidf {qTF2df' 

This expression is valid for eFi(2)f3, kFi(2)d, qTFi(2)d ^ 1 and is universal in the sense 
that is does not depend on the particular forms of the energy dispersion relations, 
transport time dependence on momentum or wave function structure. We obtain the 
familiar 2DEG T^ and d~'^ behaviour for the drag conductivity, proving that it is indeed 
a much more general result. If the metallic layers are immersed in a homogeneous 
dielectric medium, with dielectric constant e,-, one should multiply ( ITSl) by e^. Now, we 
notice that in the low temperature limit the intralayer conductivity for isotropic systems 
is given by the Boltzmann result 



e2 

(yaa = -^p{eFa)v%TF, (19) 

where the factor of 1/2 comes from the fact that we are in two dimensions, the density 
of states at the Fermi energy is given by 

PieFa) = f^^, (20) 

and that the carrier density is related to the Fermi momentum in two dimensions by 



kFa=J -. 21 

V 9a 

This allow us to express the drag resistivity in terms of the carrier densities as 
= ^ C(3) /47reo\ /4vreo\ {kpTf 
^^ e^e,2H^gW-2 \ ej J [ el J uTnTd^' ^ ^ 

It is also usual to express the drag resistivity in this limit in terms of the Fermi energy, 
momentum and Thomas-Fermi screening momentum. To do this we assume a power 
law energy dispersion relation, A (k) = Cak^"-, obtaining 

p^ = ^ ^(3)7r^ {ksTf 1 

6162 l3igi/3292 eFieF2 (kpid) {kF2d) (qpFid) {qTF2d) ' 
For drag between two 2DEG, /3i(2) = 2, gi(^2) = 2 (spin degeneracy), we re-obtain 
the known formula from [13]. For graphene we obtain exactly the same result, since 
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/3i(2) = 1, (71(2) = 4 (spin and valley degeneracy). Finally, for the case where each of 
the two layers are formed by graphene bilayers, /3i(2) = 2, gi(^2) = 4 (spin and valley 
degeneracy), we have an extra factor of |. 

3.2. The case of graphene at small interlay er distance 

Now we specialize to the case where both metallic layers are formed by single layer 
graphene, SLG, and analyze the behaviour of the drag conductivity when the layer 
separation is small, kpd ^ 1. In this situation we can no longer expand the non-linear 
susceptibilities for small q and need to consider its full dependence on q. In graphene 
the wave function overlap factor is 



cX,X' 



1 / k ■ (k + q 



■' k,k+q 2 



k + q 



(24) 



with A = +, — for the conduction and valence band, respectively. Therefore, the non- 
linear susceptibility in the low temperature limit (TT3|) reads 



where the factor of 4 comes from the spin and valley degeneracies. Equation ( !25l) is in 
disagreement with the expressions obtained in |103 both for the momentum independent 
and for the linearly momentum dependent transport time cases. However, we emphasize 
that in the low temperature limit ( l25ll holds for an arbitrary transport time. Now we 
notice that for q < 2k p, the static polarizability for graphene is given by [21] 

2kF _ 2ei 

— r - -qTF^ 

irnvp c^ 



XsLG{q < 2kp, 0) = -^ = -qTF^. (26) 



Therefore the dielectric function eRpyi(g, 0) still has the form given by (iT6l) . even if we 
do not assume that q is small. Since we have kFi(2)d <^ 1, we expand to first order in d 

Vi2{q) 



|t/i2(g,i 



eRPA(g,0) 

2 \ 2 r 

e^ ^ 



(27) 



2{2qTFiqTF2 + q{q + qrFi + qTF2)) 



{q + qxFi + qTF2Y (q + qxFi + qTF2f 

and the drag conductivity becomes 

an = ^y«^^^%^ [^^'^ (^^1' ^i^2) - rfX(^) {kF^, kF2)] , (28) 

where ag = e^ /{AneQVpK) is the fine structure constant of graphene and we have defined 
the functions 



i"«(*«.M= r , .. '"" — ^ji-^ji-4- (29) 



(g + (qXFl + qTF2)Y V ^^F1 V ^^F2 



I^'HkFi,k 



F2 



^ 2{2qTFiqTF2 + q{q + qrFi + qTF2)) L q^ L q^ 



(g + <?TF1 + <?TFl)^ V ^ 4A;fiV ^k^F2 
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Figure 1. Comparison between computed drag resistivity for graphene in vacuum 
using the full expression from equations (155)) and ([M]) . and several asymptotic limits 
as a function of the interlayer spacing. The curve Low T was computed setting y — 
in eRpA and expanding to first order in y the functions $ in p4p. The curve Low T 
+ Small d was computed using ([30]) and Low T + Large d using ((22|) . We can see 
that the agreement between the full expression and the low temperature one is good 
for most of the range, but that the small discrepancy is enhanced as d increases. 



where K 
becomes 



2m.m{kpi,kp2), and for graphene qte 



Aaqkp. 



The drag resistivity 



PD 



h 2H^ {kBTY 



-a. 



[l'^'^{kFi,kF2)-dI^'^{kFl,kF2 



(30) 



e2 3 eFieF2 

Identical result for d = has recently been derived in an independent work [12]. For 
the case where both layers are at the same carrier density kFi = kF2 = kF the functions 
X''°^ and X^^' simplify considerably and we obtain, 



I^''^{kF,kF) = l2a,--+{l 



48«3 log 1 



1 



4a, 



(31) 



+ 44(l-8ao)ao + 



l + Aa, 



+ 32ag {Ual - l) log f 



1 + 



1 



4a, 



Therefore, for drag between two SLG layers as the interlayer separation is increased, the 
behaviour of the drag resistivity changes from a linear dependence on d for grpd ^ 1 to 
a d~'^ dependence at large grpd ^ 1. In figure [1], we can see that the small separation 
expression (15U]) is only reliable for grpd ^ 0.2 and the large separation expression (12^ 
for qTpd ^ 20. If the graphene layers are immersed in a homogeneous dielectric with 
constant e^ this would correspond to kFd < 0.026^. and kFd > 2er, respectively. We 
therefore see, that these limits are not easy to be achieved experimentally. 



4. Drag in graphene: general formula and asymptotic limit 

As we have just seen, the dependence on momentum of the transport time is irrelevant 
in the low temperature and high density limits, but it will be important in general 



On Coulomb drag in double layer systems 10 

however. The dependence of the transport time on momentum in graphene depends 
on the dominant scattering mechanism. Both for strong short-range impurities 
(resonant scatterers) and Coulomb impurities the transport time depends hnearly on 
the momentum, r^^ = r° |fc|, where r° is a constant with units of length x time |l5j. 
The linear dependence of the transport time is assumed in [8] and we make the same 
assumption here. In this case, the non-linear susceptibility for graphene becomes: 

A,A' 

X ((A - A') ^ - X'q) S (eg^, - e^,^.,, + h.) . (32) 

We will follow the steps of [8] and assume that both layers are with high electron doping, 
so that the existence of the valence band can be ignored. In this case, we take only the 
A, A' = +, + contribution into account. Taking the non-linear susceptibilities at zero 
temperature, the drag resistivity can be written as 

PD = -^ — —jpY^alJ^{kFi,kF2,d). (33) 

The function J-" is defined as 



ffe»^-8.„..«E/iS./^:^.(' 



rco POO 7 ^-2dy/kFikF2X 

T{kFi,kF2,d)= dxx^ -^— -=^=Y^ : r^ 

Jo Jo smh^(y^£E^^IiyJ2\^RPAix,y)\ 

$1(X^)$2(X^ .„.. 

where the functions ^a{x, y) are introduced in Appendix B and x = q/ y/kpikp^ and y = 



q/ivrVkp^^F^). This is exactly the same expression derived in [8j if one notices that the 
function e{x,y) used there is related to eRPA{x,y) by e{x,y) = x^ {x^ — y"^) ejiPA{x,y). 
So far no approximation has been made in the sense that no asymptotic limit has 
been considered. A comparison of the different asymptotic behaviours computed in the 
previous section with the exact result, f p3|) . is given in figure [1] 

By the general arguments given in the previous section, the drag resistivity in the 
low temperature, high density, large separation limit should behave as d~^. However, in 
[8] and in this same limit a dependence of d~^ was obtained. In [10], the d~^ result is 
attributed to scaling of the vertex function v^^t^^ used in [8] with g^, while it should 
scale with q for a constant group velocity and a transport time linearly dependent on 
the momentum, implying that the carrier group velocity used in [8] depends linearly on 
momentum. We clarify that in [Bj, as in this paper, the group velocity used is constant 
and the transport time depends linearly in momentum, so that the vertex function 
depends linearly in momentum, Vj: ^tj: ^ = XvpT^k. However, as we have argued in 
the previous section, the momentum dependence of the vertex function is irrelevant in 
this limit. The incorrect d~^ dependence was obtained in [8] due to an error at the 
end of the asymptotic calculation: for small x and y with y < x the function ^{x,y) 
was taken to behave as ^{x,y) ~ -, when it actually behaves as $(x, y) ~ -^. This 
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changes the integration kernel obtained in [S] in the asymptotic hmit from q^ sinh~^ (qd) 
to g^ sinh~^ (gd), which changes the dependence from d~^ to the correct behaviour of 
d~'^. Therefore, for ei?i(2)/3, kpK^^d ^ 1, the exact ( 133|) and ( l34l) give 

^^ _ ^ C(3)vr^ {ksTf 1 ^ 

e2 2^ eirieF2 (/cfiC?) {kF2d) (qxpid) {qTF2d) ' 



in agreement with the general result discussed in the previous section. Expressing (j35j) 
in terms of the carrier density we get 

h C(3) 1 {ksTf 

^"^ e'2^^al^y^nfnfnTd^' ^ ^ 

as obtained before, but here we have started from the general expression for the drag, 
that is, fl33|) . Note that, in this limit the drag resistivity decreases as Og is increased. 
We can understand this as follows: as ag increases the screening becomes more effective 
making the momentum transfer between layers less effective. 

5. Comparison with experiments 

In the experimental setup of |3j we have two graphene layers, which we denote by t (top) 
and b (bottom). Between the two graphene sheets we have a layer of AI2O3, thickness 
d = dt = 7 nm (in [8] it was used a value of 14 nm). The bottom layer is on top of 
a db = 280 nm thickness Si02. Finally, these layers are on top of a silicon wafer. The 
carrier density in the graphene sheets is controlled using a back gate voltage between 
the silicon wafer and the bottom layer. For the relation between the gate voltage and 
the carrier densities the reader is referred to [HI IB]. Given the carrier densities, one 
can compute the drag resistivity using (!33|) and (!34|) . For the fine structure constant 
of graphene we will use the accepted value of ag = 2. To compute the drag, we must 
determine the form of the bare Coulomb interactions, Vab{q) (C]), with a,b = t,b. In a 
naive treatment of the Coulomb interactions one could write [8]: 

(naive) (^air + CAZ2O3 (naive) (naive) (^Al203 ~r £Si02 /n^\ 
^tt - 2 ^^^ ~ ^^'203 (^bb - 2 ^ ' 

with e^ir, ^AhOi and €5^02 the dielectric constants of air, AI2O3 and Si02, respectively. 
Using the values of e^ir = 1, ^Si02 =3.8 and e^ZaOs = 5.6 [22] we computed the drag 
resistivity curves with this model for the interactions. In figure [2](a), we can see the 
comparison between the experimental results from [3j and this calculation. Clearly 
the agreement is not very good. Although (1371) is correct is some limits, to obtain 
the exact form of the Coulomb interactions in layered dielectric medium one must 



solve Poisson's equation (see Appendix A ). For a 3 layered dielectric structure, for the 
effective dielectric constants eab{q) in ([7]) we would have to use the functions t^bi^) 



defined in [Appendix A[ with ei = e^jr = 1, £2 = ^Aho-i = 5.6 and €3 = esi02 = 3-8. 
In figure [2](b), we can see the comparison between the computed resistivity using this 
interactions and the experimental results. Although this is a more rigorous treatment 
for the interactions than the naive one, (1371) . the computed drag curves deviate even 
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Figure 2. Comparison between computed and measured drag resistivity. The dotted 
marks are the experimental results from [3]. The solid lines are the computed curves. 
Curves with the same colour correspond to the same temperature, (a) The computed 
curves were calculated using the naive interactions ( t^^™'^ ) . (b) The computed curves 
were calculated using the Coulomb interactions obtained by solving Poisson's equation 
in a 3 layered dielectric depicted in the inset I e^j^ (g) ) . (c) The computed curves were 
calculated using the Coulomb interactions obtained by solving Poisson's equation in 
a 4 layered dielectric as depicted in the inset ( e^^j, (g) 1 . The value used for the fine 
structure of graphene, a^, was ctg = 2 in all plots. No attempt to fit the data was 
made. 
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further from the experimental resuhs. One could proceed as in [S] and use ag as a 
fitting parameter. We found, however, that it would be necessary a value of the order 
ttg ~ 5 and ag ~ 16, for the interactions with ej^"" and e)jj, , respectively, to fit the 
data; both values are unrealistic. Taking into account that the Si02 layer is finite when 
solving Poisson's equation has virtually no effect in the computed drag resistivity. This 
situation is puzzling at least. However, more attention must be paid to how the devices 
from [3] are constructed. The AI2O3 dielectric layer is deposited in two steps, (i) First 
a 2 nm aluminium layer is deposited on top of the bottom graphene layer, being later 
oxidized. We will refer to this as AlO^; layer, (ii) On top of the oxidized aluminium 
layer, a 5 nm AI2O3 layer is directly deposited through atomic layer deposition. We will 
refer to this as ALD layer. Completed the procedure a 7 nm AI2O3 is obtained. The 
reason the AI2O3 is not directly deposited on top of graphene is due to graphene being 
chemically inert. In [22], the dielectric constant of the ALD alumina was determined 
by studying how the capacitance of devices similar to the ones from [3] scales with the 
tickness of the ALD layer keeping the AlO^ layer thickness fixed, obtaining the already 
referred value of eAhOs = 5.6. The AlO^; layer contributed with a finite capacitance, 
which corresponds to a dielectric constant of eAio^ = 2.7. The reason why eAio^ differs 
from e^iaOs is not clear but it is most likely due to interface effects between the graphene 
layer and the dielectric or between the ALD and the AlO^; layers. In this case, it is not 
clear if we can attribute a bulk value for eAio^- Nevertheless, we still considered this 
situation. For this situation of a 4 layered dielectric we use the functions ej^^ (q) from 



Appendix A , now with 62 = ^aio^ =2.7 and w = 2 nm. The comparison between the 
drag curves computed this way and the experimental data can be seen in figure [2](c). 
We can see that there is a better approximation with the experimental results, but this 
should be regarded with caution. 

6. Conclusions 

In this paper we have showed that in the limit of low temperature, high density, large 
interlayer distance and strong screening, the drag resistivity should always behave 
with T^ and d~'^. This result was obtained for general dispersion relation, momentum 
dependence of the transport time and electronic wave function structure. It is therefore 
a more general result than previous ones. Central to this fact is the general expression 
derived for the non-linear susceptibility, (fT3l) . in the low temperature limit to lowest 
order in frequency. Being general, this result also applies to graphene, and should close 
the ongoing debate regarding the behaviour of drag for this system in the aforementioned 
limit. We also derived an asymptotic expression for drag in graphene in the limit of 
low temperature, high density and small layer separation. Finally, we compared the 
available experimental data on drag in graphene from [3] with our theoretical model 
and found out that if one insists in the conventional value for graphene fine structure 
constant, a^ = 2, a quantitative understanding of the data is still lacking. This is most 
likely due to the complexity of the dielectric substrate used and one hopes that the 
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theoretical model will be quantitatively more successful in devices built using simpler 
substrates, such as boron nitride |23] . Another possibility to account for this discrepancy 
could be spacial inhomogeneities in the chemical potential of the graphene layers, which 
are not taken into account with the present formalism. Although, it is not clear if these 
should be important in the experimental range investigated, spacial inhomogeneities in 
graphene's carrier density will for sure play an important role in the electron-hole puddle 
regime and should be object of future work. 
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Appendix A. Electron-electron Coulomb interaction 

To obtain the correct form of the electron-electron interaction in a dielectric material 
one must solve Poisson's equation, — V(eV0) = Pfree/^o, where pfree is the free charge 
density. We are interested in situations where the dielectric constant is a piecewise 
constant function of the z coordinate, with discontinuities at positions Zi. Since, we 
want the interaction matrix element between states with well defined momentum in 
the X, y directions and well defined position in the z direction, it is useful to introduce 
4>{q,z) = f d'^x(f){x,z)e~^^'^, where q = {qx,qy) and x = {x,y). For a point charge 
located at {x,z) = (0,Zp), pjree = e5(0)(5(-2p), we obtain the following equation 

- ^ (<^)^'^(^' ^)) + ^''^(^' ^) = ^J(^p)- (A-1) 

The potential (j){q^z) is continuous everywhere, and the the function t{z) S^ is 
continuous except at the position of the point charge. Therefore we have the boundary 
conditions 

Solving equation f lA.l|) together with the boundary conditions (1A.2|) and imposing that 
0(g, z) decays at infinity, we can obtain the potential created by the point charge and 
from that the bare electron-electron Coulomb interaction. Let us consider that the 
metallic plates are located at z = and z = d. The electron-electron interaction can be 
cast in the form given in ([7]) with a,b = t, b, where t refers to the top layer and b refers 
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to the bottom layer. For a 3 layered dielectric: 

solving Poisson's equation gives 

e'^'^les + 62) (ei + eg) - (€3 - e2)(ei - €2] 



15 



e^^Mz) 



, z > d 

, d > z > 0, 

, > z 



(A.3) 






2(1 + e29'i) 62 - 2 (1 - e'^1'^) £3 
e^g'^(ei + e2)(e2 + £3) - (d - e2)(e3 - £2) 

4e2gde2 
e^i'^iei + €2) (€3 + €2) - (ei - 63) (eg - £2) 



(A.4) 



2(l + e29<i)e2-2(l-e29'i)ei 
This result was previously given in [9] and [21]. If we consider a 4 layered dielectric, 



e(^)(.) 
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, z > d 
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(A.5) 



we obtain, 



<:\:'(i) 






,2g(d+to) 



>2 + e'2) (63 + e',) + e''^ (62 - e',) (e', - e^ 



+6^"- (e2 - 6'2) (e3 + e'2) + e^^- {e^ + e',) (e'2 - eg 
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(A.6) 



D(^)(g) 



,2q(d+w) 



>i + £2) (e2 + e'2) - e''"^ (ei + e2) (e2 - e'^ 



+6^""' (e2 - e'2) (eg + e'2) - e''"^ {e, + e^) {e^ - e'^ 



-1 



where D^^\q) is defined as 



D('\q) 



^2q(d+w) 



ei + 62) (62 + 4) (4 + es) + e'^'' (ei + e2) (e2 - e'2) (e'2 - e^ 



+e^^" (ei - e2) (e2 - e'2) (e'2 + e,) + e'^^ (ei - e2) (e2 + e'2) (e'2 - e,) . (A.7) 



Appendix B. Details of the computation of the drag resistivity in graphene 

The contribution to the non-linear susceptibility of graphene ( 132|) coming only from the 
conductance band. A, A' = +, + is 



r(g,a;) = S-KVFr'qj — -2/^^^^^^^^^^,+) " ^^^ 



(2 
-2vpT^q[m.x^^ {q, 00) 



-k+q,- 



S I efc,+ - ^+<i,+ + ^ 



(B.l) 
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where X^^(Q') ^) is the contribution to the graphene polarizabihty coming only from the 
conductance band, 

At zero temperature, this can be computed analytically and the result can be written 



as 



Imx-^+(g» = ^J ' ^ $(g,c.), (B.3) 

1- (^ 



where the function $(g, w) is defined as 

$(g, w > 0) = $+(g, cu)e (—-q + 2A;i. ) 6 ( g - — 



+$-(g,a;)e ( A;^ - ^ - |A;^ - g| ) , (B.4) 



where 



2 



$^(g,a;) = ±cosh-^ ( '^" '=^"/"" ) T '^" ^^"/"" J (^ ^^^ ^^"/^" J -1. (B.5) 

The drag conductivity therefore becomes 

^^ - 2H^h^ksT I ^^^X sinh^ (/3na;/2) ^ _ ^^^ ' ^ ^^^ 

The intralayer graphene conductivity of layer a at low temperature is given by Caa = 
e'^vpUaT^/h = e^eFTF/ij^h^) where kpaT^ = Tpa is the transport time at the Fermi 
level. We write the interlayer interaction as f/12 = e^ ex\i{—qd)/{2eRPA{(i,oj)ei2{q)eoq). 
Introducing the adimensional quantities, x = q/ \/kpikp2 and y = uj/{yp^/kpikK2)^ the 
drag resistivity becomes 



2^62 kpT ^^7o """""" io -,inh^ (., fivFVkFikF2 

where ag is the fine structure constant of graphene. 
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